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ABSTRACT 

The optical depth is widely used in Galactic microlensing studies as a 
means to determine the density of MACHOs, since in theory it depends only 
upon their spatial distribution and is therefore less model-dependent than 
other microlensing observables. However, since the measured optical depth 
is restricted to that of events with measurable timescales, inferences of total 
MACHO density are dependent on the assumed timescale distribution. Using 
the 2-year LMC results of the MACHO collaboration, we illustrate this point 
by showing how, for an assumed isothermal halo, the inferred MACHO halo 
fraction as determined from optical depth estimates depends upon MACHO 
mass. 

The analysis highlights the following conclusions: (1) The MACHO density 
inferred from optical depth measurements depends on the assumed MACHO 
mass function for a given Galactic distribution function; (2) without extra 
information on the MACHO mass function, such as can be obtained from the 
rate-timescale distribution, optical depth measurements can provide a lower 
limit but not an upper limit on the MACHO density for a given Galactic 
distribution function; (3) a comparison between the inferred total optical depth 
of different Galactic models, or of different components in a multi-component 
Galactic model, requires knowledge of the underlying timescale distribution for 
each model or component. 

For our assumed Galactic halo model we find excellent agreement between 
our model-dependent lower-limit halo fraction of / > 0.3 (84% confidence) or 
/ > 0.15 (97.5% confidence), and the MACHO collaboration's own model- 
independent optical depth lower-limit estimates. MACHO's optical depth upper 
limits are consistent with the minimum value of the upper limit derived for our 
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assumed model. 

Subject headings: Dark matter — gravitational lensing — Galaxy: halo - 
Galaxy: stellar content — Galaxy: structure 
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1. Introduction 

Several microlensing experiments are undertaking searches to detect compact halo 
dark matter (MACHOs), as well as other low- luminosity stellar populations, following the 
suggestion of Paczyhski (|1986|) . One of the principal quantities which characterises the 
microlensing properties of a particular MACHO population is the optical depth r. This 
quantity determines the average number of microlensing events in progress at any instant 
in time per background source star. In its simplest form it is given by 

t = / — -p(x) dx, (1) 
jo m 

where x is the distance along the line of sight between the observer and MACHO, L is the 
observer-source distance, m is the MACHO mass, p is the MACHO mass density at x and 



AGmx(L-x) 

R* = \l (2) 



is the Einstein radius. Equation ([!]) is valid if all sources are at distance L, otherwise one 
needs to further integrate equation ([I]) over the spatial distribution of sources ( |Kiraga & 
Paczyhski 19941 ). 



So far, theoretical predictions have been compared to measures of the optical depth 
obtained from observations towards the LMC and Galactic bulge ( |Alcock et al. 1997a| ; 



Renault et al. 19971 ; |Udalski et al. 1994j ; |Alard et al. 1995| ; |Alcock et al. 1997b| ), with 



further tentative comparisons also starting to emerge from observations towards SMC 
( |Alcock et al. 1997cj ) and from pixel experiments directed towards M31 ([Ansari et al. 1997 



|Crotts fc Tomaney 1997| ). For non-pixel based experiments, observational determinations 



of the optical depth are based on the model-independent estimate 



1 7T ^ s t. 



e.g. |Udalski et al. 199~^ ; |Alcock et al. 1997a| ; |Alcock et al. 1997E| ), where t e j (i = 1 . . . N ohs ) 



are the measured event timescales (we define t c = 2R e /Vr, with Vr the MACHO velocity 
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across the observer-source line of sight), £ is the efficiency with which timescales t e are 
detected and E is the "effective exposure" ; that is the average observation time per source 
star multiplied by the total number of stars observed. Equation (^) basically measures the 
fraction of the total observing time for which microlensing events are in progress. Errors 
are typically determined by a boot-strap method in which the range in optical depth is 
estimated from random timescale realisations generated from the observed t e ^ (e.g. |Alcock 



et al. 1997b| ; [Alcock et al. 1997a| ). Han & Gould ( |1995[ ) have shown that error estimates 



based on naive Poisson statistics can significantly underestimate the true error. 

Alcock et al. (1997a) point out that equation (^) is not a measure of the total optical 
depth, but only of the optical depth of events which fall within a particular range of 
timescales (those for which £ > 0). They stress that an estimate of the total optical depth 
requires one to input a timescale distribution. 

In this study we define the concept of observable optical depth and use it, together 
with the MACHO collaboration's 2-year LMC results, to place limits on the total optical 
depth, and hence MACHO halo fraction, for the "standard" isothermal halo model analysed 
by MACHO. The analysis highlights the dependency of the results on MACHO mass for 
a given Galactic model. Whilst optical depth measurements allow one to determine a 
lower limit on the density of MACHOs for a given Galactic distribution function, no upper 
limit can be obtained without using extra information, such as can be obtained from a 
rate-timescale analysis. 



2. Optical depth-timescale distribution 

Because of the timescale dependence of equation (|3|), its theoretical analogue is not 
equation ([!]), since this is an implicit integral over all event timescales. Instead, one can use 
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the following expression: 



7T r dr . 

r = - t e —dt e , 4 

4 Jt e (£>o) dt e 



where T is the event rate. Equation (H) is almost a re-statement, in somewhat expanded 
form, of the relation r = (7r/4)(t e )r, where (t e ) is the average event duration. The one 
difference is the restriction to timescales t e (£ > 0). 

Equation (f|) points to an expression for the differential contribution to the optical 
depth from events of duration t e : 

dr 7r dr 

dT e = 4 te d%- (5) 

Evidently, just as the rate-timescale distribution dr/dt e depends upon the spatial, velocity 
and MACHO mass distributions, so too must the optical depth-timescale distribution 
dr/dte- Its integral, if performed over a restricted range of timescales, is similarly model 
dependent, unlike the expression in equation (TI]). Hence, to evaluate r within a certain 
range of event durations t e one must specify both the full Galactic distribution function 
and the MACHO mass function. 

As an example, we employ the cored isothermal halo model originally analysed by 
Griest ( |1991| ), and denoted model S in the MACHO collaboration's halo analyses. For the 
case of a discrete mass function and stationary line of sight, the rate-timescale distribution 
towards the LMC for this model is 

dr _ 2V 2 p Q/ 2 2 f L f3(m } x) 2 exp[-(3(m,x)} 

a + -n-oJ / T~h o„d „„„;,„„„; , „2 , nl\ ax \°> 



dt e m Jo (x 2 — 2xR Q cos b cos I + a 2 + Rq 

(c.f. Griest |1991| ), where V c = 220 km s" 1 is the halo velocity normalisation 



Po = 0.0079 Mq pc~ 3 is the local halo density, a = 5 kpc is the halo core radius, 

Ro = 8.5 kpc is the Sun's Galactocentric distance, (I = 280°, b = —33°, L = 50 kpc) is the 

LMC position in Galactic coordinates and (3 = (2R e /V c t e ) 2 . From equations @ and (|6|), 
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the optical depth-timescale distribution becomes 

dr _ nV 2 p t c ^ 2 | ^ 2 r L p(m,x) 2 exp{-P(m,x)] ^ 
dt e 2m Jo (x 2 — 2xi?o cos 6 cos Z + a 2 + Rq) 

Both distributions are plotted in Figure (0) for the LMC direction, assuming a MACHO 
mass m — 1 Mq. To ease comparison, the distributions are both normalised such that their 
peak values are unity. It is evident from the figure just how much more sensitive the optical 
depth is than the rate to longer duration events. 



3. Observable and observed optical depths 

Having defined the optical depth-timescale distribution we can now compare theoretical 
prediction with observation in either of two ways. The first method would be to directly 
apply equations (|) and (£|) to microlensing data. Instead, we proceed by defining the 
observable optical depth 

Toble^ / £(t e )^dt e . (8) 
Jt e (£>o) at c 

This quantity represents the total optical depth which is potentially observable to a 
microlensing experiment with detection efficiency £. Since it already incorporates the 
detection efficiency, it should be compared not to equation (|3[) but to the directly observed 
optical depth in the absence of efficiency corrections: 

1 7T ^ s 
i=l 

In the limit of low-number statistics, which is presently the case for searches towards the 
LMC and SMC, comparison of the quantities r bi e and r b s d should provide a more stable 
estimate than can be obtained from equations (|]) and @), since it is less dependent on the 
efficiency estimates for particular timescales. 
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Fig. 1. — The optical depth-timescale and rate-timescale distributions towards the LMC for 
a standard isothermal halo and a lens mass m — 1 Mq. The line of sight is assumed to 
be stationary with respect to the Galactic rest frame. The peak of both distributions is 
normalised to unity for ease of comparison. 
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An evaluation of the quantities r bi e and r b s d allows a straightforward estimate of the 
halo fraction / = T b B d/7"obie f° r our adopted halo model. This fraction represents the total 
MACHO fraction, not just the fraction within some timescale range. The first 2 years of 
MACHO observations towards the LMC have uncovered 8 microlensing candidates ( [Alcock 



et al. 1997a|) , including a likely binary event (Pennett et al. 1996| ). MACHO analyses both 



this sample and a six-event sub-sample which excludes the binary event, since it seems 
likely that the event originates from within the LMC, and also excludes one of the weaker of 
the other candidates to preserve the overall average event duration. We restrict our analysis 
to this 6-event sub-sample. The timescales of these 6 events yield an observed optical depth 
r obsd = 5.7 x 10~ 8 for an effective exposure E = 1.82 x 10 7 star-years. Note that this is 
smaller than the value quoted by MACHO since it does not compensate for the effect of 
efficiencies. Instead, it represents the optical depth actually measured by the experiment. 

As mentioned previously, the potentially observable optical depth r bi e is sensitive not 
only to the halo model but to the assumed MACHO mass function. We calculate inferred 
halo fractions for discrete MACHO mass functions with masses ranging from 0.01 — 10 Mq. 
This range includes the 0.1 — 1 Mq range favoured by the MACHO collaboration's 
maximum-likelihood analysis of the rate-timescale distribution ( Alcock et al. 1997a ). The 



MACHO 2-year detection efficiencies are incorporated as required in equation (BJ). 

The solid line in Figure (|]) shows the preferred halo fraction / as a function of assumed 
MACHO mass m. The errors on / are determined by Monte-Carlo simulation and are 
discussed in the following section. The most obvious point to note is the mass dependency 
of /. Whilst the preferred value for / is quite stable between ~ 0.1 — 1 Mq at / ~ 0.5, it 
becomes arbitrarily large for small and large MACHO masses, with a minimum preferred 
value of / = 0.47 occurring for m = 0.23 Mq. This variation is directly related to the 
timescale correspondence between the efficiency and optical depth distributions. For very 
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Fig. 2. — Halo fraction determinations resulting from the 2-year MACHO LMC 6-event 
sample. A discrete MACHO mass function is assumed. The solid line denotes the preferred 
value as a function of m. Also shown are 68% and 95% confidence-level regions bounded 
by dashed and dotted lines, respectively. The horizontal dashed line at / = 1 demarcates 
a full MACHO halo, whilst the cross denotes the MACHO collaboration's 2-D maximum 
likelihood solution based on the rate-timescale distribution (m = 0.41 Mq, / = 0.51). 
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high- or low-mass MACHOs the peak in the optical depth-timescale distribution occurs at 
relatively long and short timescales, respectively, where the efficiency is low. Hence, r bi e is 
very small for these cases, and so / oc r~^ le becomes large for a given r b s d- Hence, for a 
particular Galactic distribution function, optical depth measurements alone can specify / 
only as a function of m. 

4. Error estimates 

To obtain error estimates on our inferred halo fraction / we use a Monte-Carlo 
procedure in which 10000 microlensing "experiments" are conducted for each assumed 
MACHO mass m and halo fraction /. The errors are obtained from the resulting 
distribution of "observed" optical depths. 

The calculation procedure is as follows: for each assumed mass m the efficiency- 
corrected rate-timescale distribution £dY /dt e and optical depth-timescale distribution 
£dr/dt e are calculated for a full MACHO halo (/ = 1), together with their integrals r o bi e 
and T bie- Comparison of r bi e with r b s d = 5.7 x 10~ 8 , as computed from the MACHO 2-year 
LMC results, gives the preferred estimate of f(m) discussed in the previous section. For a 
range of assumed /, the expected number of detectable events for the model is computed as 
N exp (m) = fETobic. This number is used to generate a Poisson realisation iV b Si j for mass 
m, fraction / and experiment j. Event durations £ e ,i (* = 1 . . • -Wobsj) are generated from 
the distribution SdT/dt e . The optical depth observed in experiment j, Tj, is then computed 
from equation (^|). Repeating this process for each experiment results in the distribution 
P Tj (f,m), of which some fraction F Tj>Tohsd (f,m) give Tj in excess of the measured value of 
5.7 x 10~ 8 . For a given (/, m), the quantity (Tj), the average over all experiments, is found 
to be typically in agreement to within 0.5% of the quantity /r bi e as obtained by direct 
integration over £dr/dt e - 
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For an assumed m, the values of / for which F T . >Tobsd is 0.16 and 0.84 (0.025 and 0.975) 
bound a 68% (95%) confidence region in /. These confidence intervals are bracketed by 
the dashed (dotted) lines in Figure (Q). It is clear that these intervals are larger for larger 
lens masses m. This is because a measured r b s d implies a fixed total time J2i t e ,i fo r the 
summed event durations. The number of events required to produce this summed duration 
is inevitably larger for low-mass MACHOs than for high-mass MACHOs, due to the smaller 
Einstein radius of low-mass MACHOs [c.f. equation (0)]. As a result, N exp is required to be 
large for low-mass MACHOs so Poisson fluctuations about A eX p are small. The converse is 
true for very massive MACHOs, which give rise to inherently longer durations. Since A exp 
for these objects is typically smaller, larger Poisson fluctuations can arise, which in turn 
produce larger errors in /. 

The actual number of detected events, A obs , is of course a known quantity but is not 
directly used in the optical depth analysis. This is because the analysis is concerned only 
with the sum of the event timescales J2i t e ,i not the number of events, and this is essentially 
why the analysis places no constraint on the MACHO mass m. Rate-timescale analyses 
use the observed number of events together with the probability of observing each event 
duration (for an assumed underlying timescale distribution) to constrain both the MACHO 
density and mass. 

Figure @ shows that, using optical depth measurements alone, one can obtain firm 
mass-independent lower limits on / for the assumed halo model. The measured value for 
7"obsd = 5.7 x 10 -8 implies lower limits of / > 0.3 (84% confidence) and / > 0.15 (97.5% 
confidence) for the model adopted here. Both of these estimates are in excellent agreement 
with the MACHO collaboration's model-independent boot-strap determination ( |Alcock et 



al. 1997a]) , and thus support MACHO's assertion that the lower limits it derives ought to 



be robust. Figure @ also highlights the fact that upper limits cannot be placed without 
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extra information on the MACHO mass. The MACHO collaboration's boot-strap analysis 
gives upper limits which are consistent with the minimum derived upper-limit values in 
Figure (@). These minimum values are / < 0.7 (84% confidence) and / < 0.9 (97.5% 
confidence). Note that these upper limits apply only to specific MACHO masses, and that 
the mass at which the upper limit minimises is a function of the required confidence level. 
The comparison of upper-limit estimates confirms the MACHO collaboration's concerns 
that its model-independent boot-strap method may tend to underestimate upper limits, 
since it does not take account of the possible contribution of events with durations exceeding 
the experiment lifetime. 

Plotted in Figure (Q) is the MACHO collaboration's preferred value for / and m 
based on a 2-D maximum likelihood analysis of the rate-timescale distribution (|Alcock 
et al. 1997a| ). The correspondence between this estimate (/ = 0.51, m = 0.41 Mq) and 
the optical depth f(m) constraint for the same mass is reassuring and shows that optical 
depth measures can also serve as important consistency checks on rate-timescale analyses. 
One would not have the right to expect such good agreement in cases where the assumed 
model is a poor approximation of the actual Galactic distribution function, or where the 
underlying timescale distribution is poorly sampled by the experiment. For m = 0.41 Mq 
we obtain / = 0.48]to3i, where the quoted errors bound a 95% confidence interval. This 
implies a total halo optical depth r = 2.241^45 x 10~ 7 , which is to be compared to the 
MACHO collaboration's own model-independent estimate of r| 00 = 2.06i?Jg x 10~ 7 for 
events with durations between 2 and 200 days ( |Alcock et al. 1997a| ). Comparison of the two 
values shows that the MACHO collaboration's estimate is slightly lower than ours, which is 
to be expected since our estimate is not restricted to a certain timescale range. However, 
the MACHO estimate is not much lower because, for the assumed Galactic model, the 
timescales produced by 0.41 Mq lenses typically fall within the MACHO efficiency range 
and are thus relatively well sampled. 
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5. Discussion 

The optical depth is a familiar concept in gravitational microlensing studies and is 
widely used as a means to determine relatively model-independent constraints on the 
density of MACHOs in our Galaxy In this study we emphasise that a comparison between 
the observed optical depth and the predicted total optical depth is a model-dependent 
procedure and that failure to take proper account of this may give rise to misleading results. 

We have shown that even if one specifies a distribution function for the MACHO 
population, optical depth constraints on their density will still be a function of the MACHO 
mass. Whilst optical depth measurements can be used to place firm mass-independent lower 
limits on the MACHO density they cannot be used to place upper limits. The rate-timescale 
distribution, however, can provide both lower and upper limits on the MACHO density 
because of its ability to simultaneously constrain both the MACHO density and typical 
mass. 

In cases where more than one Galactic component contributes significantly to the 
observed lensing rate, as is believed to be the case for searches directed towards the Galactic 
bulge, a proper calculation of the relative contribution of each component to the observed 
optical depth requires one to assume both distribution functions and mass functions for 
each component. Such assumptions are implicit in calculations which simply compare 
observed and theoretical optical depths without regard to timescale distributions. This 
is also relevant to LMC searches if the LMC or an intervening structure is contributing 
significantly to the microlensing statistics ( |Sahu 199~4] ; |Zhao 199~T| ; ^aritsky fc Lin 1997[ ). 



Similarly, if one wishes to make a comparison between the halo fractions of different Galactic 
models, as inferred from optical depth measurements, on must specify full distribution 
functions and mass functions for each model. 

In summary, optical depth measurements can provide firm lower limits on the density 
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in MACHOs and can serve as useful consistency checks for rate-timescale analyses. By 
themselves, optical depth measurements cannot provide upper limits on the MACHO 
density. Comparisons between the preferred MACHO density of different Galactic 
components, or of different Galactic models, are meaningful only if one assumes both a full 
distribution function and MACHO mass function for each component or model. 

The author wishes to thank Will Sutherland for providing the MACHO 2-year LMC 
efficiencies and David Bennett for pointing out errors in an earlier version of this work. 
This research is supported by an EU Marie Curie TMR postdoctoral fellowship. 
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